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ABSTRACT 



A Bayesian likelihood analysis of the stellar density profile shapes of several dark 
matter dominated Galactic dwarf spheroidal galaxies (dSphs) has been completed. 
This analysis indicates that the inner power law, p oc r -7 , of each of these profiles 
satisfies 7 < 1.21 at greater than 99% confidence with a slight preference for 7 ~ 0.7. 
Speculating that this bound on the stellar profiles may be related to the inner power 
law of the underlying dark matter distribution reveals a statistical discrepancy with 
the values favoured by recent N-body simulations of cold dark matter (CDM), j C dm ~ 
1.4- 1.5. 

Key words: methods statistical - galaxies: dwarf- galaxies: haloes - galaxies: struc- 
ture - dark matter. 



1 INTRODUCTION 

The nature of the dark matter remains one of the prime 
areas of astrophysical and cosmological research. There has 
been some controversy recently concerning cold collisionless 
dark matter (CDM), arguably the simplest solution to dark 
matter problem. Assuming that satellite dark haloes can 
be associated with dwarf spheroidal galaxies (dSphs) , CDM 
appears to overpredict by two orders of magnitude the num- 
ber of these satellit e haloes surrounding the Galactic halo 
jMoore et al. 1999af) . Recent N-body simulations of CDM 
also indicate that dark matter haloes have cuspy cores on 
subgalacti c scales, with an i nner expone nt Icdm ~ 1.4 — 1.5 
iFukushige fc Makino 19971: IMoore et al. 1999b : 

Jing fc Suto 200(1 ISubramanian et al. 200* : 



Bullock et al. 2 001). 



up fr om the value found in earlier 
- 1.0 jNavarro. Frenk. fc White 19961 



simulations of j C di 
hereafter NFW). 

Recent detailed, high resolution observations 
of the rotational motion of LSB galaxies indicate 
their core density profi les are su bstantially flatter 
than CDM prediction s llde Blok et al. 20011 : see also 
iBorriello fc Salucci 200 lfl. (For earlier discussion and de- 
bate, see for exampl e Ivan den Bosch fc Swaters (200CT) and 
IMoore et al. (1999bD ). 

Since the discrepancy between CDM theory and obser- 
vation is greatest on small scales, it is natural to seek to 
sharpen it further by examining the smallest, nearest galax- 
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ies, namely dwarf spheroidal galaxies. Galactic dSphs have 
long been recognised as excellent probes of the nature of 
the dark matter. Besides being the smallest known galax- 
ies they also possess the highest known M/L ratios. Ob- 
servationally, their proximity means that individual stars 
can be resol ved, so that beam sme aring effects of the type 
discussed by Ivan den Bosch fc Swaters (20001) for LSBs are 
not an issue. Unfortunately, whilst velocity measurements 
have established the presence of large amounts of dark 
matter, they have hitherto not been sufficiently accurate 
to determine the radial profile of the dark matter density 
JArmandroff. Prvor fc Olszewski 1997) . 

H owever, the stellar density profi le is rather well estab- 
lished lllrwin fc Hatzidimitriou 1995T) and enough data ex- 
ists for a detailed comparison with the theoretically pre- 
dicted dark matter profiles. This comparison is only mean- 
ingful if the stars trace the mass, a very strong assumption, 
which is hard to justify. Nevertheless, it seems a reasonable 
first step to fit dSph stellar profiles to the parametric form 
used in theoretical studies of dark matter profiles, whilst ac- 
knowledging that any conclusion on the dSph dark matter 
profile must await more accurate velocity measurements. As 
both theory and observations develop, we anticipate that 
a more complete likelihood analysis, of the type performed 
here but including both density and velocity data, should 
become possible. 

One should also point out that at least for some star for- 
mation histories, comparing the stellar and dark matter pro- 
files is not without meaning. If the bulk of the stars formed 
before the dSph dark matter haloes, they would thereafter 
cluster under gravity just like dark matter particles. If, on 
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the other hand, stars formed after the dark matter haloes, 
from baryons which were heated when the Universe was 
reionised at z > 5, or alternatively if the stars in dSphs 
were accreted from other parts of the galaxy, there would 
not necessarily be any simple relation between the stellar, 
and dark matter density profiles. 

CDM theory over-predicts both the number density, 
and the velocity dispersion of haloes with the masses of 
dwarf spheroidals jBode. Qstriker. fc Turok 20001) . If CDM 
is to be viable, there must exist some mechanism (such as 
stellar feedback) to suppress star formation in most, but 
not all, low mass haloes. The surviving low mass galaxies 
would then be those in which stars formed earliest, and those 
which possessed the highest central density. These selection 
effects would tend to favor haloes with steeper-than-typical 
inner profiles and higher central densities. Comparing a typ- 
ical CDM halo profile against observed dark matter profiles 
would actually be conservative in this respect. 

In a similar way we may also speculate on how com- 
patible other theories of dark matter are with observa- 
tions of dark matter dominated dwarf spheroidals. For ex- 
ample, it has been suggested that a warm dark matter 
(WDM) scenario may better match the observed num- 
ber of satellite galaxie s predicted for the local group 
JAvila-Reese et al. 20011) and so a mechanism to suppress 
star formation may be less necessary. Additionally, stars may 
be a more viable tracer of the dark matter in WDM theo- 
ries as the recollapse of small WDM haloes is expected to 
occur considerably later than in the CDM case. However, de- 
spite these advantages, theories of WDM may predict haloes 
cuspy enough to be cause for concern as in the CDM case 
feke et al. 200 ill 



2 METHODOLOGY 

We have studied the radial distribution of stellar mat- 
ter in the dark matter dominated Galactic dSphs 
using the framework of Bayesian likelihood analysis 
jGreeorv fc Loredo 1992^ . 

The stellar density profile was assumed to take the 
smoothl y interpolating double power law form of equation 
(1) fe.g.. lKravtsov et al. 199gh . 
Po 



p(r) 



(r/r )T[0.5 + 0.5(r/r ) a ](f 3 ~i)/< : 



(1) 



Using this form yields p oc r -7 for radial distances 
smaller than ro, and p oc r _/3 for radial distances larger 
than ro. The parameter a controls the scale over which the 
power law shifts from 7 to j3. Roughly speaking, this change 
occurs between radii e~ 1 ' a 'ra and e 1 ' Q To. This density pro- 
file accommodates a wide r ange of use f ul models, includi ng a 
modified King-like profile <King 19621 : iRood et al. 19721) for 
(a, 13, 7) = (2, 3, 0), the NFW profile for (a, /?, 7) = (1, 3, 1), 
and the profile suggested by recent N-body CDM simula- 
tions (a, f3, 7) = (1, 3, 1.5). 

Projecting onto the celestial sphere gives the effective 
two dimensional stellar surface density. 



E(r,) = 2 



r 2 , + z 2 ) dz 



(2) 



Currently, the best observational determination of stel- 
lar density as a function of radius in the dSphs of interest 



is the work of llrwin fc Hatzidimitriou f 19951) . and can be 
summarised in this analysis by a vector E (, s which contains 
the mean density at a discrete set of radii and an associated 
covariance matrix C h s . 

For each model (a, (5, 7,ro,po) we have calculated the 
conditional relative likelihood £(a,[3, 7,ro,po|Eobs) -af 
statistic in practice - to quantify the compatibility of a par- 
ticular model with observations. 

The chief advantage of the Bayesian formalism is that 
it allows one to assess the relative likelihood of a particular 
subclass of the models by integrating £ over the relevant 
equiclassification surfaces in the model space. To assess the 
relative likelihood of different values of 7 we integrate C 
over all models with the same 7 but different (a,/3, ro,po) 
and obtain the marginalized distribution £ 7 (7|E & S ). This 
process allows for a fair decoupling of our state of knowledge 
of 7 from (a,/3,ro, po) without, for example, ignoring the 
data in the outer regions of the dSph in an ad hoc manner. 

Care must be taken when choosing the ranges of param- 
eter values to ensure that a fair sampling of the parameter 
space is made into all regions that can contribute signifi- 
cantly to the marginalized distribution. The reason for this 
is the edges of the parameter distribution impose the triv- 
ial prior that the parameters of candidate models must lie 
within the range of parameters chosen. For the purposes of 
this study we have found that a largely unbiased range of 
parameters is: a G [0.5,12]; (3 G [2.5,17]; 7 G [-13,3]; 
r G R; and p G R. 

In practice C is an extremely rapidly decreasing func- 
tion of (ro, po) about its maximum, and so for computational 
purposes the marginalization over these variables has been 
approximated by a maximization. We find that maximiza- 
tion over (ro, po) yields indistinguishable results to marginal- 
ization provided that the characteristic turnover radius, ro, 
is not permitted to be unreasonably large (ie. larger than 
the edge of the galaxy). Marginalization over (a, j3) proceeds 
straightforwardly to yield £ 7 . 

The three dSphs we have focused our analysis on, Ca- 
rina, Draco and Ursa Minor have mean ellipticities, e = 
1 — b/a, of 0.33, 0.29, and 0.56 respectively and the avail- 
able surface density data is averaged over elliptical bins. To 
make a conservative analysis in spherical theory we made a 
direct comparison of circularly rebinned surface density data 
with the theoretical profiles, using the properly reweighted 
error distribution. In general, we find that even the rough 
approximation of comparing the elliptically binned density 
profile with spherical theory gives similar limits on 7 as the 
gross dependence of C on the shape parameters is insensi- 
tive to the exact data binning scheme. This is expected as 
spherically smoothing an ellipsoidal distribution has the ef- 
fect of smoothing the radial distribution over scales of ~ ea 
without significantly altering the gross radial dependence - 
especially in the nearly flat inner regions important for this 
analysis. 



3 DISCUSSION 

The marginalized likelihood distribution C-, for the three 
dSphs in question is shown in Figure 1. Although the peak 
of the distribution fluctuates slightly from galaxy to galaxy 
the upper bound on 7 is consistent between all three. 
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Figure 1. The likelihood distribution of 7, the exponent of the 
inner power law, for Carina (dot-dashed), Draco (dotted), Ursa 
Minor (dashed), and for the three galaxies taken together as an 
ensemble (solid). 



While the distributions in 7 have tails stretching well 
into the negative region of parameter space this is not un- 
expected. Negative values of 7 are simply models that have 
a peak in the density at a scale of ~ ro, a plausible result 
given the adopted uncertainties. Projection effects coupled 
with the freedom to renormalize and rescale the density pro- 
file have the result of extending the tail relatively far into 
the negative 7 region. Despite the ambiguity between mod- 
els in the negative 7 regime, it is clear from the Figure 1. 
that models with 7 > 1.21 are strongly excluded by this 
analysis. 

The best-fitting density profiles are shown in Figure 2. 
Note that while the profiles are nearly always good fits for 
the inner and central regions there sometimes appears to be 
an excess of density for the outer few data points - related 
to the scheme a dopted for background subtraction of stars 
(Irwin & Hatzidimitriou 1995). The results with respect to 
the parameter 7 do not change if these points are included 
or excluded from the analysis. We have adopted the strategy 
of including all available data. 

Table 1 summarises the limits this analysis places on 7 
and the other shape parameters for each galaxy on its own, 
and the galaxies considered together as an ensemble. It also 
presents the limits with a tophat prior of 1 ^ a ^ 2 imposed 
on the likelihoods - the range of parameter space that the 
NFW and recent N-body CDM profiles are contained in. The 
effect of this prior is to shift the peak of the distribution and 
the upper bound on 7 to even lower values than when the 
full range of a is used. 

As measurements of the radial profile of the velocity 
dispersion in these galaxies increase in sensitivity a joint 
analysis fitting both the stellar dispersion and stellar density 
profiles predicted by a given DM profile should provide a 
much better test of the CDM and similar paradigms. 



4 CONCLUSION 

We have completed a Bayesian likelihood analysis of the 
stellar density profile shapes of the Galactic dwarf spheroidal 




Figure 2. Shown are the maximum likelihood models (a, [3, 7) = 
(5.5,4,0.7) for Carina (Top), (a,/?, 7) = (9.5,3.5,0.6) for Draco 
(Middle), and (a, j3, 7) = (3.5,4.0,0.5) for Ursa Minor (Bottom). 
Models with less sharp turnovers of a ss 1 — 2 can also be found 
that yield reasonable fits. 



Table 1. Summary of shape parameter limits. 



Uniform Priors 



Galaxy 



7(99%)* y P e a k P (2<r) a (2a) 



Carina 


<1.09 


0.80 


4.0i 2 j 


<5.2 


Draco 


<1.12 


0.60 


,5+1.9 


<7.3 


Ursa Minor 


<1.21 


0.80 


4-0± ;l 


<2.7 


Ensemble 


<0.98 


0.70 


4.0i°;« 


- 




1 ^a sg2 


Tophat 


Prior 




Carina 


<0.76 


0.00 


5.0+t; 


- 


Draco 


<0.76 


0.00 


5.0T?° 4 


- 


Ursa Minor 


<0.95 


0.00 


B.0151 


- 


Ensemble 


<0.52 


0.00 


™±li 


- 



*These values are based on the one sided Gaussian positive devi- 
ation from the peak. Calculating with the full two sided distribu- 
tion yields an even stronger upper confidence bound on 7. 

galaxies Draco, Ursa Minor, and Carina. With uniform prior 
assumptions for the parameters (a, /3, 7, rrj,po) we conclude 
that, for each galaxy, 7 < 1.21 at greater than 99% confi- 
dence. 

As discussed in the introduction, except in the special 
situation where the bulk of the stars formed before the dark 
matter halo, jnu would not be expected to be the same as 
7 stars. Nevertheless the comparison we made quantifies the 
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the discrepancy between the observed inner stellar profile 
of dSphs and simulated dark matter haloes, a discrepancy 
which CDM theory or one of its cousins must ultimately 
explain. 
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